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1. Introduction 

For some time now it has been standard practice in e + e~ reactions to compare event- 
shape distributions with resummed perturbative predictions (see for instance [1]). 
The resummation is necessary because in the two-jet limit (small values of the shape 
variable) the presence of large logarithms spoils the convergence of the fixed-order 
calculations. Such resummed analyses have led to valuable information about the 
strong coupling constant and also about non-perturbative effects [2]. 

At HERA, similar studies of DIS event-shape distributions are being carried 
out by both collaborations [3,4], but as yet no perturbative resummed calculations 
exist for comparison between data and theory. Here we present the results of such a 
calculation. 

So far resummed predictions for processes involving initial state hadrons exist 
for observables such as DIS and Drell-Yan cross-sections in the semi-inclusive limits 
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where x for DIS and r for Drell-Yan are close to one [5-7]. In addition there are 
predictions for jet-multiplicities in DIS [8] and the W transverse momentum in Drell- 
Yan production [9] which entail somewhat similar considerations to those that are 
encountered here. 

In DIS, event shapes are defined in the current hemisphere of the Breit frame 
to reduce contamination from remnant fragmentation, which is beyond perturbative 
control [10]. The distribution of partons in the current hemisphere is analogous to 
that in a single hemisphere in e + e~, enabling us to adopt the usual e + e~ methods 
for the resummation. But it turns out that event-shape variables also depend on 
emissions in the remnant hemisphere (typically through recoil effects) — a proper 
resummed treatment of the space-like branching of the incoming parton is therefore 
required. 

The fact that only one, predefined, hemisphere is used in DIS event shapes leads 
to a number of subtleties compared to e + e~ event shapes. Event shapes measure 
properties of the energy-momentum flow in the current hemisphere. They are always 
dimensionless, so it is necessary to normalise this measure of energy-momentum flow. 
Two choices are possible: a normalisation to Q/2, or to E, the energy in the current 
hemisphere (whereas in e + e~ the total energy E is equal to Q). Furthermore, one 
often requires the choice of an axis along which to project the momentum flow (for 
example for the thrust, the broadening). In e + e~ there is one logical choice of axis, 
namely the thrust axis (that which maximises the thrust). In DIS, two axes arise 
naturally: the photon axis and the thrust axis (which is given by the direction of the 
sum of all momenta in the current hemisphere). 

In this article we shall consider two variants of the thrust. Both are measured 
with respect to the photon axis: one, Tq, is normalised to Q/2 and the other, T E , 
is normalised to E. In the limit of a 1+1 jet event both tend to T = 1. Since we 
are particularly interested in this limit, it will be convenient to define r = 1 — T and 
consider the limit of small r. While we are generally interested in the differential 
distribution it will turn out to be more convenient to look at the integrated 
event-shape cross section 



from which the distribution can be straightforwardly obtained by differentiation. 

For small r one finds that the shape cross-section contains terms a™ ln m r, m < 
2n coming from multiple soft and collinear radiation. The lnr factors compensate 
the smallness of a s , hence the need arises for resummation, which leads to a result 
schematically of the form S(r) ~ exp(— G^ats In 2 r + ...), with G±2 some number 
which will be calculated. From the point of view of classifying the terms that we shall 
calculate, it turns out to be most convenient to consider InS. This contains terms 
a" ln m r with m < n+1. Terms with m = n + 1 are referred to as leading logarithmic 




(1.1) 
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(LL), while those with m = n are known as next-to- leading logarithmic (NLL). We 
shall resum both these sets of terms, and neglect as subleading the remaining terms, 
m < n. 

A word of caution is needed about te, and in general about all event-shape 
measures which are normalised to the energy in the current hemisphere, E. At O (a s ) 
there are configurations in which the current hemisphere is empty. At higher orders 
it can be filled by soft-gluon radiation from the partons in the remnant hemisphere. 
If the event-shape is normalised to E then it can take on finite (significantly different 
from and/or 1) values due to the presence of soft-radiation; in other words it will 
be infrared unsafe starting from O (a 2 ). This has been known for some time. The 
standard experimental practice is to exclude events in which the energy in the current 
hemisphere is too low. For example, HI set the threshold E > £u m = Q/10 [3]. 
Though formally this eliminates any infrared unsafety, one is still left with potentially 
large logs a s (a s lnQ/2£ lim ) n . Thus we recommend that £\ im be chosen somewhat 
larger, say £\ im = Q/A. Though this will slightly reduce the statistics, it ought to 
improve the convergence of the perturbation series. 

The structure of this article is as follows. In section 2 we discuss the kinematics 
and define the thrust measures which will be studied. In section 3 we perform the 
resummations. Section 4 gives the expressions to be used for the practical calculation 
of resummed distributions. The results are expanded to second order and compared 
to the predictions from the fixed order Monte Carlo programs DISASTER++ [11] 
and DISENT [12]. This enables us to cast some light on the differences between 
DISASTER++ and DISENT, observed in [13]. Finally a matching prescription is 
proposed for combining the resummed and fixed order predictions. This is followed 
by the conclusions. 

The appendices provide details about the leading order calculation of the thrust 
(appendix A), the techniques used for the inverse Mellin transforms (appendix B) 
and the integrations needed for a result correct to NLL (appendix C). 

2. Thrust definitions and kinematics 

It is convenient to write the momenta ki of radiated partons (gluons and/or 
quark- ant iquark pairs) in terms of Sudakov (light-cone) variables as (Fig. 1) 

ki = aiP + foP' + k u , Oifc = k 2 ti /Q 2 (2.1) 

where P and P' are light-like vectors along the incoming parton and current direc- 
tions, respectively, in the Breit frame of reference: 

P = xP p , P' = xP p + q, 2(P ■ P') = -q 2 = Q 2 (2.2) 

where P p is the incoming proton momentum and x = Q 2 /2(P p ■ q) is the Bjorken 
variable. Thus in the Breit frame we can write P = ^(5(1,0,0,-1) and P' = 
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tj<5(1, 0, 0, 1), taking the current direction as the z-axis. Particles in the current 
hemisphere Hq have $ > ai while those in the proton remnant hemisphere 7Yr have 
Oil > Pi. 

The momenta p and p' of the initial- and final-state 
quarks can also be resolved along the Sudakov vectors P 
and P' . From momentum conservation, 



p= (l + a' + ^2a^j P, 

p> = a 'P+(l-J2 Pi) P '-J2 k * > ( 2 - 3 ) 

wnere Figure 1: Kinematics. 

a/(l- J>) = fchtf /Q 2 . (2-4) 

The initial quark is assumed to be collinear with the proton direction and is therefore 
aligned along P (neglecting possible 'intrinsic' transverse momentum). 




2.1 Thrust normalised to Q/2, T Q 

Consider first the thrust Tq, defined by 
current hemisphere normalised to Q/2: 




the sum of longitudinal momenta in the 



where the index a run over all partons and the the index i runs only over emitted 
partons. We have assumed that the outgoing quark momentum p' lies in the current 
hemisphere, because the consideration of the non-soft, non-collinear emissions that 
would result in this not being true is beyond our concern in the present case. We 
have 

T Q = l-T Q = Y t Pi- " ^) + a ' 

H C 

= p* + ai + a ' = min i a i, p^ + a ' ■ ( 2 - 6 ) 

Wr He 

Where not explicitly stated, the sums extend over all emitted partons, irrespective 
of which hemisphere they are in. Here the contribution from those partons that are 
not directly in the current hemisphere arises through the recoil effect that they have 
on the current quark. 
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2.2 Thrust normalised to E, T E 

Another way to define the current jet thrust is to normalise it to the total energy in 
the current hemisphere, instead of Q/2. This gives 

T E = T Q /(l-e), (2.7) 

where 

e= 1 - | E E - ■ ( 2 - 8 ) 

^ a&ic 

The quantity e, called the energy deficit in the current hemisphere, is itself an inter- 
esting shape variable [14]. It is given by 

He Wr He 



Hence, we have, 

te — 1 — Te 



3. Resummation 

We shall consider a reduced integrated cross section, where by reduced we mean 
that it will be the contribution to F 2 (rather than to da/dxdQ 2 ) which comes from 
configurations with 1 — T < r. There will be no need to consider contributions to 
Fl since they are suppressed by powers of r for small r. 

We write the following expression for the (reduced) cross section for the emission 
of m gluons in the remnant hemisphere and n in the current hemisphere, where all 
the gluons have k t <C Q, 

e 2 a (O 2 ) 1 FT r ^M)Cf d 1 A [ Q2 dkl 3 a s {kl)C F 



with 



f 1 + z 2 

pR,i = / dzi Q(Qai - k u ), (3.2a) 

J I - Zi 

pc, = J ^ (i + (i - m em - k t] ). (3.2b) 



N is the moment variable (conjugate to Bjorken-x) and 



an = ^ — ^- . (3.3) 

Z\ . . . Zi 
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We have introduced a factorisation scale Qq, and accordingly in consider only 
emissions above that scale. Qq is chosen so as to be much smaller than any other 
scale in the problem. 

For the event shapes considered here, the independent gluon emission pattern 
written above is enough to reproduce the behaviour of the full QCD matrix elements 
up to the required (NLL) accuracy, with the following restrictions: a s must be eval- 
uated to two loops in the CMW (or Physical) scheme [7, 15]; we must also put in the 
virtual corrections and take into account all possible branchings of the incoming leg, 
not just the emission of gluons from a quark; but these can both be done trivially 
later on. It is not necessary to take into account the splitting of emitted partons, 
other than what is already implicitly included through the running of the coupling. 

The fundamental difference between the expression given above and the corre- 
sponding one for e + e~ comes from the presence of the quark distribution. Hard 
collinear emissions change the momentum fraction of the incoming parton, and this 
must be accounted for, through a (multiple) convolution of some function (to be 
determined) with the parton distributions. We take moments (N) with respect to x 
in order convert that convolution into a product. These moments enter only for the 
emissions in H R — emissions in He are essentially identical to those in a hemisphere 
of e + e~. 

The use of moments does not however remove all difficulties. In particular in 
eq. (3.3) we are still left with a product of z^'s. This matters both in the function 
in eq. (3.2a) and in the contribution to the thrust, which for say tq will be = 
ku/(Q 2a i) ( c f- ec l- (2-6)). The solution relies on the fact that the gluons which 
contribute to the thrust are either the gluon with the highest transverse momentum 
or soft gluons. Considering the gluons as ordered in transverse momentum, the 
gluon with highest transverse momentum is number 1 and ol\ = (1 —Z\)jz\. One can 
equally well consider the gluons to be ordered in angle, and if gluon m has the highest 
transverse momentum, then gluons with % < m must have 1 — Zi <C 1 and one can 
write a m ~ (1 — z m )/z m . For soft gluons we can approximate the splitting function 
by 2/(1 — Zi), and then replace (1 — Zi)j{z\ . . . 2j_i) =>• (1 — Q). The integration 
measure remains 2dQ/(l — Q), and ctj just becomes (1 — 0)/0- So we can in general 
just replace a, by (1 — z^jzi (where Zi should sometimes be understood to mean Q). 

3.1 Resummation of tq 

To write a resummed expression we take the Mellin transform of our G function for 

e(-£ft-E«<-«')=/^(nH (H' )' (3 ' 4) 

This almost factorises the expression for the function into pieces which each depend 
on a single emission. There remains the less trivial term e~ ua ': from (2.4) we know 
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that ol is approximately (using 1 — ^/3 ~ 1) the squared vector sum of emitted 
transverse momenta. For now it is actually sufficient to consider the sum of squared 
momenta. For soft and collinear gluons k^/Q 2 = aiifa is negligible compared to 
mm(ai, Pi). Only for collinear but hard gluons is k 2 JQ 2 comparable to min(a i , / 9j) 
because one of the ctj, fa is O (1). In the current hemisphere, since fa < 1, the 
contribution from the k 2 /Q 2 can be neglected: the region where it is comparable to 
«j is single-logarithmic and its inclusion modifies r by a numerical factor, i. e. changes 
a s lnr to say a s ln2r, a difference which is NNLL and so negligible. The difference 
is similarly NNLL in the remnant hemisphere, however the fact that can be much 
greater than 1 (as large as 1/x), and hence k^/Q 2 ^> fa, means that one should at 
least keep track of it. Therefore in the remnant hemisphere contribution we retain 
a contribution to r of the form k 2 /Q 2 (keeping track of it more accurately would 
involve using the squared vector sum — but for tq the error is down by one more 
logarithm). Thus we replace 

We then sum over m and n in (3.1) to give us the following exponentiated form 
for the integrated thrust distribution 

Siv(r) = e 2 q q N (Q 2 ) J ^- e -^M-^M . (3.6) 

The current hemisphere radiator, TZc is 

KcM = _ / f^L £ dp l±<^l &m _ kt) {e - m _ lh (37) 

where the —1 accounts for virtual corrections (as in [16]). The remnant hemisphere 
radiator is a little trickier. The attentive reader will have noticed that eq. (3.1) 
has a leading factor of qN(Ql), whereas eq. (3.6) has qN(Q 2 )- Thus 1Zr(v) should 
contain a (//-independent) piece to take into account the change in scale of the quark 
distribution, 

1n Qn(Q 2 ) f Q2 dk 2 a s (k 2 ) 
Qn^Qo) JqI kf 2tt 

(where 7 gg is the standard quark anomalous dimension) as well as the part arising 
from the sum over m in eq. (3.1). For now we are still working in a framework 
involving only gluon emission from a quark, hence the use of only the 7 gg part of the 
anomalous dimension matrix to change the scale of the quark distribution. 
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Therefore we have 

" Q2 dk 2 a s (k 2 )C F 



2 k 2 t 2tt 



^ j 1 dz l -±^-Q{Qa - k t ) (yr e ->>W+W<P) _ i) 



-2^2), (3.9) 

where as before we have introduced a term —1 to account for virtual corrections. 
Next, we replace a — (1 — z)/z. In the function we can neglect the z in the 
denominator (the 6 function is relevant only for small 1 — z). So we obtain 



dk[a s (kf)C F ( f ^l + z 2 
IT 



I dz 

7 i- 



7 99 (iV) \ 



(3.10) 



We note that in the limit of small z the combination (f3+k 2 /Q 2 ) just reduces to k 2 /Q 2 , 
i.e. a limit on r is just equivalent to a limit on the emitted transverse momentum. 
Then we write 

7«W = CV jT 1 - 1)9(1 " * ~ kt/Q) + O (|) , (3.11) 

to obtain 

~ Q2 dk 2 a s (k 2 )C F f 1 



(3.12) 

The result no longer depends on our original choice of factorisation scale, Qq, and 
thus we drop the lower limit on the kt integral. Our next step is to make the 
approximation, valid to our accuracy [16], 

e -Q^fe) _ ij„ _ e _ g2( *~^ (3.i3) 

with P = z/e 7E . The fc t integral can then be divided into two pieces (according to 
whether the limit on 1 — z in the 6-function of (3.13) is above or below 1) 

J Q 2 9 -i kf 2tt J 1-z 



Jq2 P -2 k t 2n j\-vk\ic? 1 - ' 
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where in the second term we have made the approximation 1 — z <C 1. We express 
the z integral in the first term as 



A 1 + z \n n , Z- 1 -^. 1 + z 2 lqq (N) Q 3 

dz- -(z -1)+/ dz- ~ = —^ — + 2 In---, (3.15) 

1 — z J 1 — z Of k t I 



where we have neglected terms of order kt/Q. The z integration in the second term 
of (3.14) gives 2\n(Pk t /Q). 

The current hemisphere radiator can be evaluated in a similar manner, but with- 
out the complications from iV dependence. The sum of the radiators is 



k 2 2tt \ C F ' ""k t 

f Q2p ^ dk 2 a s (k 2 )C F A , , , , 
+ / TT oj F 41n(z7A; t /g). (3.16) 



We see that the quark anomalous dimension j qq (N) appears in our answer. It arises 
from the restriction of transverse momenta of hard emissions. A consideration of all 
branchings leads to r y qq (N) being replaced by the full anomalous dimension matrix. 
This is a standard result in a situation where a restriction is imposed on the final-state 
transverse momentum [9] . The inclusion of all possible splittings leads one, however, 
to consider also soft gluons emitted from a gluon. Fortunately they can be ignored in 
the present situation because it can be shown that for such a gluon to contribute to 
the thrust it would have to be emitted at a larger angle than the last hard emission, 
and thus, by coherence [17] it is emitted with the colour charge of a quark. To see 
that the angle really is larger, let us consider two gluons contributing equally to the 
thrust (in the remnant hemisphere), i.e. with similar (3 values. The angle of the gluon 
with respect to the incoming quark is proportional to (3/k t . Gluons with transverse 
momenta smaller than that of the last hard emission and contributing equally to the 
thrust are bound therefore to have a larger angle. 
Thus our final answer in v space is 



^ e 2 q N (Q 2 /v) exp 



de t a s {e t )c F f Aln Q_ 3 



Q2p -i k 2 2ir \~~'k, 



- r 1 §^ iH > kt/Q) 



(3.17) 



where the integral involving the anomalous dimension matrix has been absorbed into 
a rescaling of the quark structure function. For brevity we rewrite the contents of 
the square bracket as —2R(u). 



9 



The inverse Mellin transform is performed using the techniques of [18], as out- 
lined in appendix B and gives 



^n(t q ) ~ e 2 q q N (T Q Q 2 ) 



T(1 + 2R') 

1 



77 exp 



/e 7E 
-2R — 



T(l + 2R>) 



77 exp 



-2R 



1 



- 2 lE R' 



(3.18) 
(3.19) 



where R 1 = R'(1/tq). In the scale of the quark distribution, numerical factors 
such as e 7E have been thrown away since they correspond to NNLL effects. Explicit 
expressions for R and R' are given in appendix C. 



3.2 Resummation of t e 

Starting from the expression for te, (2.10), we have the O-function 







r 



He 



^2 hi 



IQ 2 



(3.20) 



We have taken the contribution to a' only from transverse momenta in the remnant 
hemisphere, because all transverse momenta in the current hemisphere contribute 
much less (except for emissions with (5 ~ 1 and a ~ r, but this region is not 
enhanced by any logs and so can be neglected). 
We introduce 



and write the 6 function as 




(3.22) 



Analogously to what was done in the previous subsection, we write the integrated 
cross section as 

= e 2 q g N (Q 2 ) J ^e™/ 2 J ^ e ^e*V^M . (3.23) 

Since the contribution from partons in the current hemisphere is the same as in the 
tq case (except for the factor of 2 included in the definition of the Mellin transform), 
IZc is the same as in that case, i.e. 1Zc(v) = Rip). 
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Arguments analogous to those used for tq, eqs. (3.8) to (3.12), lead us to the 
following result for the remnant-hemisphere radiator TZ R (b), 



n R (b) ~ - 



c F r Q2 d% 



2tt 



Tlk? 



a s (k t ) dz z 



1 + z 2 
1-z 



Q(l-z-k t /Q) e 



-ib.kt 



- 1 



(3.24) 



£f$* W (2^ +2h g-») (1 - W (3, 5) 



From [16] this integral can be evaluated to the required accuracy by replacing 

e 7E 6 



'i-j (bkt))=>e(bkt-i), 



b = 



(3.26) 



Thus we have 

TZ R (b) = Ru(b*Q>) + r ^^M) 7w (JV), (3.27) 

J 1/52 fc t 27T 

where Ru(v) is defined in appendix C. We now perform the p-integration in (3.23) 
to get 

'1/6 2 &t 



-Ru(b 2 Q 2 ) - / -rf 7«W - Rip) 

(3.28) 



± ] I ^_L p ~b 2 Q 2 /4u I u _ 



~ e " qN \Mi J ^2 



R'r 



-R v {bl)-R{u) 



(3.29) 



where we have expanded Ru(b 2 ) = Ru(pfy +ln& 2 /b^R^ipQ) + (a s ) and substituted 
7 9g with the full anomalous dimension matrix and used it to change the scale of the 
quark distribution. The integration over b yields 



^n(v) ^ e 2 q q N v .,, 



b 2 Q 2 ^ R ' u 
4i> 



r(i-^) e -^( 5 o Q2 )-^ p ), 



(3.30) 



For simplicity we then choose b\Q 2 = Au, or equivalently b^Q 2 = e jE u. 

The inverse Mellin transform with respect to v is similar to the that for tq, and 
we obtain as our final answer (remembering that v was conjugate to r/2) 



^n(t e ) ~ e 2 q N (T E Q 2 



e 2 q qN(r E Q 2 



ra-R' 



U) 



T(l + R! + R! v ) 

, r(i-i^) 
r(i + R' + r'jj) 



cxp 



cxp 



Te 
1 

Te 



1 \ 



(3.31) 
(3.32) 



( 7 E + hi2)#-(2 7E + ln2)i^ 



As before, R' = R'(l/r). 
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4. The final result 

It is possible to go one step further in the resummation, namely the determination of 
constant terms at O (a s ). These are terms independent of r which essentially multi- 
ply the resummed answer. Their function is to compensate for mismatches between 
the independent gluon emission pattern (as applied to a single gluon) and the actual 
emission pattern, usually in corners of phase space, regions not enhanced by any logs. 
The corners that can arise are the region of hard, non-collinear emission, and, for 
the thrust, the region of collinear, non-soft emissions with k 2 ~ rQ 2 . There are also 
other origins for such terms, such as certain approximations made in the expression 
for the thrust, and contributions from virtual factorisation scheme corrections. 

The contribution to F 2 from events with 1 — T < r is therefore given by the 
following expression: 



and we have written the non-x dependent leading and subleading logarithms as 
Lgi(a s (3oL) and g2(a s j3oL) respectively, whose forms follow from eqns. (3.19) and 
(3.32) and the expressions for R and Ru in appendix C. Specifically, for tq we have 
(where it should be noted that the Ri and R 2 are defined with different arguments 
from R) 




(4.1) 



where L = In 1/r and the convolutions are defined as 




(4.2) 






(4.3b) 



and for t e , 

#i(a s /3 £) 
g 2 (a s (3 L) 



R 2 (a s p L) 



Rui(&sPoL), 
Ru2(a s (3 L) 



(4.4a) 
(4.4b) 

r(i - r' v ) 



( lE + \n2)(R' + R' u )- lE R' u + \n 



r(i + R' + r' v ) ' 



It is useful to have their order by order expansions, 




(4.5a) 



n 




(4.5b) 



n 
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where the O (a s ) and O (a 2 ) coefficients are given in table 1. 





G\2 




C23 


C22 




-2C F 


3C F 


27rf3 Gi2 






-\c F 


(3-31n2)CV 


16-7T a r~t 
— P0<J"12 


l^ 2 C F + C; x l 22 \n2)C A C F 
+(l + l\n2)n f C F 



Table 1: Coefficients in the resummation. 



The expressions for the C±s are given, in the DIS and ms factorisation schemes, 
in appendix A. These constant terms have been associated with a scale Q 2 , both 
in a s and in the parton distributions. Strictly such an identification of the scale 
is not unique: indeed C\ can come from pieces with a variety of scales, however 
the error made in neglecting this fact amounts to terms a™ +l L n and thus is beyond 
our accuracy. A further point worth mentioning relates to the cutoff on the visible 
energy, £i im . For both tq and t f , C\ is independent of £n m , as long as events with 
less than E < Su m in the current hemisphere are given a thrust of 0: this is because a 
modification of £\\ m just redistributes thrust configurations between some finite value 
of t and t = 1; but the small r region is not affected. (For certain other variables, 
such as the thrust measured with respect to the thrust axis and normalised to E, 
this would not be the case). 

4.1 Comparison with NLO computations 

A valuable test of the resummation is to expand our answers to second order in a s 
and compare the result with fixed order predictions from programs such as DISAS- 
TER++ [11] and DISENT [12]. Though agreement does not guarantee the correct- 
ness of the resummation, it is nevertheless a non-trivial check. 

For the comparison we simply look at the second order contribution to 

1 dZ N (Q 2 ,r)^ 



a 0:N (Q 2 ) rflni ' 

where a is just the O (a") part of the structure function moment F 2j n(Q 2 )- Moment 
space with respect to x is used to reduce the various convolutions (for the Ci's and 
the change of scale of the parton distribution) to products. At order a 2 , the terms 
that we control go as a 2 L 3 , a 2 L 2 and a 2 L and their coefficients can be determined by 
expanding (4.1). Therefore the difference between our expanded resummed answer 
and the exact second order coefficient is a term which for large L is at most a 
constant. Figure 2 shows the expanded resummed and exact second order results 
(from the fixed order Monte Carlo program DISASTER+- 1-) in the quark and gluon 
channels. 
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The shape of the resummed and exact results are clearly similar and the difference 
between them does seem to tend to a constant at small r. This is visible more clearly 
in figure 3 which shows the difference between the two results. Also shown is the 
difference between the current version (0.1) of DISENT and the resummed answer. 
In the gluon channel the difference is clearly inconsistent with a constant: it seems 
to go as a 2 L 2 . In the quark channel there is also evidence of similar problems, with 
a significance of a few (between 3 and 5) standard deviations. Results are shown for 
te in figure 4, where the same pattern of successes and problems is seen. 

The origin of the discrepancy in DISENT has yet to be elucidated. However it 
has been possible to identify that the disagreement lies in pieces proportional to Cp 
(quark channel) and CaTr (gluon channel). Furthermore the fact that, in the gluon 
channel, the disagreement is at the level of a a 2 L 2 term in the distribution (a 2 L 3 
in the integrated distribution) suggests that the problem lies with a term involving 
one soft and two collinear divergences. In the gluon channel there should be no 
term a 2 L 2 with colour factor CaTr (due to angular ordering, soft gluon radiation 
in the region relevant to the thrust is associated with the colour factor C F because 
it is radiated at an angle larger than that of the quark emitted in the photon-gluon 
fusion). 

Above we have said that the difference between the resummed and DISAS- 
TER++ is consistent with a constant. This should be qualified: in the gluon channel, 
the largest-allowed non-constant is quite small, about ±0.2(a s /2n) 2 L for both tq and 
te- In the quark channel the statistical error on the results is much larger, because 
the largest term goes as a 2 L 3 (in the gluon channel it goes as a 2 L 2 ), and for a given 
number of events the statistical error is proportional roughly to the largest term. 
The largest allowed non-constant term here would go as ~ ±3(a s /2n) 2 L. In both 
the quark and gluon channels the uncertainty represents about 10% of the coefficient 
of the (a s /27r) 2 L term from the resummation. 

The DISASTER++ results have been obtained 1 with the equivalent of 50 days' 
time on a machine with a SPEC CFP95 [20] of about 20 (or roughly equivalently, 
500 MIPS). Thus significantly higher precision is not really feasible. (DISENT, for 
comparable errors, is an order of magnitude faster, so when the the problems with 
DISENT have been resolved a more stringent comparison will be feasible; DISENT 
also shows less dependence on the internal cutoff, required in order to avoid floating 
point errors). 

4.2 Matching 

The above form for the answer is satisfactory for very small r. However in practice 
one wishes to fit to data over a range of r. In e + e~ sophisticated methods have been 

1 Using the Condor system [19] to distribute the processing across INFN computers spread across 
the whole of Italy 
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Figure 2: For tq, a comparison of the coefficient of (a s /27r) 2 from the expanded resummed 
answer and that from DISASTER++. Shown for N = 1.4. 
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Figure 3: For tq, the difference between the second order coefficient from the Monte 
Carlo programs (DISASTER and DISENT) and that from the expanded resummed answer. 
Shown for N = 1.4. 



developed for combining the O (aj) and the resummed calculations in order to have 
an expression which can be applied over the whole range of data [1] . This procedure 
is usually called matching. It essentially involves two elements: the removal of terms 
which would otherwise lead to double counting (i.e. logarithms present in both the 
fixed-order and resummed answers) and the exponentiation of the subleading log- 
arithms (a^L which is associated with the coefficient G21 of g^(a s (3QL)) which are 
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Figure 4: For te, the difference between the second order coefficient from the Monte 
Carlo programs and that from the expanded resummed answer. Shown for N = 1.7. 



present in the fixed-order calculation. 

In DIS the situation is a little more delicate because the term a 2 L which would 
normally be associated with G21 is difficult to exponentiate — for example some of it 
might arise because the parton distributions convoluted with the CYs should really 
have been evaluated at a different scale, or because of contributions from the NLO 
splitting functions. Some of these contributions exponentiate only in moment space 
(with respect to x) and as a matrix in flavour space. Thus the simple exponentiation 
carried out in e + e~ is not appropriate here. 

However such an exponentiation is not strictly necessary. What is necessary [1] 
is to ensure that terms such as a 2 L are not independent of the resummation, but 
rather a subleading modification of the resummation; thus our final answer must 
contain not a 2 L on its own, but a\L exp(Lgi(a s PoL) + ...). 

As yet we have not yet carried out a detailed study of possible matching schemes. 
However for completeness we do include a preliminary proposal for matching, which 
satisfies the requirements discussed above and which will be examined in more depth 
elsewhere [21]. 

Since these matched expressions will be directly compared to data, it is useful 
at this point to switch to considering the cross section, differential in x and Q 2 , for 
1 — T to be less than r. For the resummed differential cross section, which we will 
represent by cPctr, we have 

, _ d 2 a R (x,Q 2 ,r) 47ra 2 (l + (l- j /) 2 ) 2 
d aR = dxdQ 2 = ~x& 2 Q ' T) • (47) 
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For brevity it will turn out to be useful to define the following matrices: 



q(x) 



( Qu(x) \ 
qa(x) 

V 9(x) J 



P(x) = 



fP„(x) ■■■P qg (x)\ 
P„(x) 



(4.8) 



and 



c , (rc) _w ( i + (i-,) 2 ) 



Cl{x) 



Q 4 



4na 2 (1 + (1 - yf) 
Q 4 2 



fe 2 J(l-x)\ 
£5(1 -x) 



V o y 

e'Ci, 9 (x) 



(4.9) 



(4.10) 



With this notation, and using eqs. (4.1) and (4.7), the resummed differential cross 
section becomes 

d 2 a R (x, Q 2 , t) = [Co <g> q(x, rQ 2 ) + a s (Q 2 )d ® q(x, Q 2 )] e L 9 i(M) +92 (a s /3 L) ^ 

(4.11) 

where we have introduced a s = a s /2n. 

The form that we propose for the matching is the following 



d 2 a(x, Q 2 , t) = d 2 a R + [a s (d 2 a% ] - d 2 a^ 

+ a 2 (d 2 a { £ - d 2 af - - d 2 a^)(L 2 G 12 + LG u j) 



(4.12) 



where d 2 a^ is the coefficient of a™ in the exact result (as determined from Monte 
Carlo programs) and d 2 a^ is the coefficient of a" in the resummed answer: 

d 2 o { £ = [C (G 12 L 2 + G U L) - C ® PL + Cj ® q(x, Q 2 ) , (4.13) 
d 2 af = [C (l/2 G 2 2 L 4 + (G 23 + G U G 12 ) L 3 + (1/2 + G 22 )L 2 ) + C <g> (-G 12 PL 3 
+ (i P <g> P - G n P - 7vP P)L 2 ) + d(G 12 L 2 + G n L)} q(x, Q 2 ) . 

(4.14) 



For these expressions to be correct, in eq. (4.11) the evolution of the parton distri- 
butions to scale rQ 2 must be carried out from scale Q 2 with just the leading log 



17 



DGLAP equations (in any case this is all that can be guaranteed by the resumma- 
tion procedure). If one wants to use the NLL DGLAP equations, one may, but then 
eq. (4.14) should be modified accordingly. 

Finally we point out that the effect of power corrections on the distribution 
should just be to shift the distribution to larger r by an amount equal to the shift 
of the mean values as calculated 2 in [14,26], as is the case in e + e~ [23,27,28]. 

5. Conclusions 

We have determined the resummed distributions for two definitions of thrust, in the 
region T — > 1. For tq the answer is almost identical to the e + e~ result, except that 
the answer is multiplied by a parton distribution, which is evaluated not at scale Q 2 
but at scale tQ 2 . This arises because placing a limit on the thrust amounts, in the 
remnant collinear region, to placing a limit on the emitted transverse momentum. 
The thrust normalised to the energy in Tic a l so has this dependence on the quark 
distribution at tQ 2 . But it differs from tq in that the coefficient of the leading double 
logarithm, G\2 is only 3/4 of that of tq, the reason being that the two normalisations 
depend differently on emissions in the remnant hemisphere. 

These resummed expressions have been expanded and compared to the predic- 
tions from the fixed-order Monte Carlo programs DISASTER++ and DISENT. We 
find reasonable agreement with DISASTER++, but a definite discrepancy compared 
to DISENT. The discrepancy is most visible in the gluon channel, but there is evi- 
dence that there might be a problem also in the quark channel. 

We have also given a 'matching' prescription for combining the resummed and 
fixed order predictions: it satisfies the basic properties of being correct both to second 
fixed order and for the leading and sub-leading logarithms, and additionally of not 
having 'large' (constant or logarithmic) pieces left over at small r. 

There remain a certain number of other DIS event-shape variables that can be 
resummed. They include the broadening, the C-parameter and the thrust measured 
with respect to the thrust axis and the jet mass. These will be examined elsewhere 
[21]. 
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A. Leading order calculations 



Here we give the explicit expressions for the functions C\ q (z) and C\ g (z) entering the 
result, eq. (4.1), for the integrated thrust distribution. Essentially, they are obtained 
by carrying out the full leading-order calculation, taking the limit of small r, and 
then subtracting the pieces with a logarithmic dependence on r. 

We give results for the F 2 part only, in accord with eq. (4.1), since the component 
of the cross section contributing to Fl is suppressed by a power of r at small r. In 
the DIS factorisation scheme, the results for tq are as follows: 



C™(z) = —Tj 



R 



1 + z 2 , 1 l + 2z-6z 2 

In - H ; r 

1-z z 2(1 -z)+ 



(z 2 + {l-z) 2 ) In --(l-Qz(l-z)) 



(A.l) 
(A.2) 



The plus prescription has its usual meaning. The equivalent results in moment space 
(with respect to z) are 



^Dis r< 
u iq,N — 



U lg,N — ~ 1 R 



N+l N 2 (N + 1) 2 2 



+ - MN) + 7 ) + W(N) 



1 1 

+ - — + 



(A.3) 
(A.4) 



_(N + 2) 2 (N + l) 2 N 2 N N + l N + 2 

where ip(N) = d/dN\nT(N). 

The analogous results for te are somewhat more complicated, because of the 
more involved nature of the phase space integrations for that variable: 



C™{z)=-C, 



1 + z 2 



1 — z z 



In - + (1 + z 2 



f\n(2(l-z))\ 1 + A\n2 + 2z - Qz 2 



1-z 



2(1-*) 



+ 



'is 



R 



-(1 + z) ln2 + - (ln 2 2-ln2) 5(1 - z) 
2 

( z * + (i _ z f) i n + (_i + 6z _ 6z 2 ) 



(A.5) 
(A.6) 



The results in moment space are 
<- 2 2 3 



^yDIS 



71 



6 , N , N+i +^(N) + ^ + 1 4>(N)+ 1 ^(N + 2)+ 1 2 + 



+ T7 + 



+ \W(N) + + 2)) + \mN)f + (m + 2) 2 )) + 
+ 2 - In 2) - In 2 ( 2 (^(N) + 7) + 1 + 



^DIS 



2 



N N + l 

._L (1 + tP( N ) + 7 ) + _A_ (3 + ^(N + 1) + 7) 



(A.7) 



N + 2 



(3 + V(iV + 2)+ 7 ) + hi2 - 



1 



2 2 
+ 



N N + l N + 2 



(A.8) 
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Additionally to go from the DIS to the ms factorisation scheme result, it suffices 
to add the O (a s ) Ms-scheme F 2 coefficient function to the above results: 



riMS _^dis | n 



ln(l -z) 

1-z 



2(l-z) + 



-(1 + z) ln(l - z) 



^±^lnz + 3 + 2z-f^- + l)S(l-z) 



r<MS _^DIS , rp 



{(l-zf + z 2 ) In 



1-z 



- 8z 2 + 8z - 1 



(A.9) 
(A.10) 



B. Inverse Mellin transforms 

In [18] the following operator technique was developed to help carry out the Mellin 
transforms. We write 



e -n(u) = e -tt(e-«a) y -a 



a=0 



The inverse Mellin transform is given by 



E(r) 



2mv 

Using the identity 



r(l + a) \t 



a=0 



-n(e-aa) a 



x a g(a) 



a=0 



= e- n ^ d ^g(a) 



a=0 



we can absorb a power of (1/r) into the argument of the radiator 

S(T) = e U ) WTa) 
Performing the logarithmic expansion of the radiator: 

-11 (xe~ da ) = -7Z(x) + K'(x)d a - \1Z"(x)d 2 a + . . . , 



a=0 



(B.l) 



fB.2) 



(B.3) 



(B.4) 



(B.5) 



we obtain that the action of the operator on a function which is regular in the origin 
reduces to substituting TZ'{x) for a, while lZ"(x) = O (a s ) and higher derivatives 
produce negligible corrections: 



r(i + ft') 



n' = n'(-) . 



(B.6) 
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C. Two-loop radiator integrals 

There are two integrals which have to be evaluated: 

Ruiy) 

and 



1 dk 2 a s (k 2 )C F f^Q 2 3 

Q 2 u 



k? 



2tt 



fc? 2tt 



(C.l) 



(C.2) 



'Q 2 u- 2 

For the results to be correct at two-loop accuracy it is necessary that a s be the CMW 
or Physical scheme coupling [7] and that the running of a s be taken into account to 
two-loop level. The CMW scheme coupling is related to the ms coupling through 

and the two-loop running of the coupling is reproduced to sufficient accuracy by the 
expression 



K = Cj 



a s (/i 2 ) = a s (Q 2 ) 
where 



1 



A a s (Q 2 )ln(l-A) 



1-A A) 



;i-a)- 



\ = a s (Q 2 )(3 \n%, (C.4) 



11C A -2n f 17 C\ - 5C A n f - 3C F n f 

12^ ' fl ~ 24^2 • 

We split the results into leading and subleading logarithms: 

Ru{v) = In v ■ R V1 (a s (Q 2 )(3 In v) + R U2 {a s (3 Q In v) 

and analogously for Rl and R = Rjj + Rl- The results are as follows: 

C F 



'2tt/3 X 



[-A-ln(l-A)] , 



R U2 (X)An(l-X) + ^ F[X + {1 - X)Hl - X)] 



4tt/5 



+ 



C F (3 1 



4n 2 [3 2 (l - A) 
_A + ln(l-A)_l ln ' 
1-A 2 v ; 



(C.5) 
(C.6) 

(C.7a) 
(C.7b) 



and 



Rli(X) 



C f 



2ir{3 \ 



A + (1 -2A)ln 



1-2A 
1-A 



XC F [(1-A)ln^-A] C F/3l 

ttL2\*) = ; TTT^T. ~ r 



4n 2 P 2 (l -A) 



2tt/? 3 



A + (2A-l)ln(l-A) 
T^A 



(C.8a) 
(C.8b) 



+ ^ ln 2 (l - 2A) - - ln 2 (l - A) + ln(l - 2A) 
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and for the sum R = Ru + Rl, 
C F 



Ri(X) 



[(1 - 2A) In (1 - 2A) - 2(1 - A) In (1 - A)] , 



^^Ind-A)^-^^-^-^ 1 - 2 ^ 



4tt/5 
2tt(3 2 



4ir 2 (3 2 



(C.9a) 
(C.9b) 



+ 



1 



ln(l - 2A) - 2 ln(l - A) + - ln 2 (l - 2A) - ln 2 (l - A) 

2 



We also need the expressions for the derivatives of the i?'s with respect to lni/ 

C F A 



R', 



u 



2%f3 1-X' 



R'r 



R' = 



2 In 



C F 
2tt/5 

C F . 1 - A 



1-A 



A 



1-2A 1-A 



In 



7ip 1-2A' 



(C.lOa) 
(C.lOb) 
(C.lOc) 



where A = a s /3o In v. A final point relates to scale changes. If in In vR\[X) we change 
the scale of a s from Q 2 to [i 2 , then R 2 gets modified as follows: 



R 2 ^R 2 + X\n^-(R'-R 1 ) , 



(C.11) 



and analogously for Ru and Rl. 
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